Role of the superexchange interaction in magnetic transition and polaron crossover 
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The Hubbard-Holstein model is studied including double-exchange interaction and superexchange 
interaction using a variational phonon basis obtained through the modified Lang-Firsov (MLF) 
transformation followed by the squeezing transformation. The kinetic energy, polaron crossover 
and magnetic transition are investigated as a function of electron-phonon (e-ph) coupling and elec- 
tron concentration for different values of antiferromagnetic superexchange interaction (J) between 
the core spins. The polaron crossover, magnetic transition and the suppression of ferromagnetic 
transition with J are discussed for the model. 
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Recently the interest in the double exchange model 
[1] has grown considerably with the discovery of very 
large negative colossal magnetoresistance (CMR) [2] and 
anomalous magnetotransport properties in doped man- 
ganites [3], namely in Ri^A^MnOa (where R=La, Pr 
and A=Ca, Sr, Ba). Ferromagnetism in these compounds 
(for 0.2-0.4) is believed to be due to the 'double- 
exchange' (DE) mechanism which operates when local 
Mn-ion spins are strongly coupled by Hund's rule with 
the spins of the itinerant electrons occupying a narrow 
band. However, the experimental results [3] in mangan- 
ites, namely the sharp change in resistivity near T c and 
the physics of CMR, cannot be explained by the DE alone 
[4]. Millis suggested the lattice polaron effects due to 
strong electron-phonon (e-ph) interaction as a necessary 
extension [5]. Roder et al. [6] also examined the com- 
bined effect of e-ph interaction and DE on T c using the 
variational wave function techniques. In fact, the con- 
tribution of the lattice polaron to carrier mobility was 
pointed out earlier by Goodenough [7], Several theoret- 
ical models have been proposed based on lattice-carrier 
coupling [6,8-11]. Many experiments [12] indicate evi- 
dence of strong lattice-electron coupling in Manganese 
Oxides [12] which shapes the properties of manganites 
crucially. Moreover, small to large polaron crossover is 
reported by many experimental groups [13,14]. There 
are models [8,9,11] which incorporate DE interaction in 
a polaronic model. Min and co-workers [8] studied the 
combined model of spin double exchange and lattice po- 
laron to investigate the effect of small to large polaron 
crossover on the magnetic and transport properties under 
the mean- field approximation scheme. However, the su- 
perexchange interaction also plays an important role in 
manganites and a study considering e-ph and superex- 
change interactions in a double exchange model is im- 
portant. 

In manganites, a competing tendency towards the lo- 
calization comes through the e-ph interaction which sta- 



bilizes a local distortion of the Oxygen octahedron sur- 
rounding each Manganese ion. DE favours the ferromag- 
netism whereas superexchange interaction (SE) antifcrro- 
magnetism; hence the DE and SE interactions compete 
with each other. The interplay between lattice-carrier 
coupling, DE and the SE interactions is the reason for 
the rich phase diagram of the manganites. So it would 
be of great interest to study the role of SE interaction( J) 
in the combined model of DE and e-ph interaction. In 
this work we will follow a method based on variational 
phonon basis [15] which is promising for Holstein [16-18] 
and related models [19] to predict the correct results up 
to intermediate range of hopping (t) (t < loq where lvq 
is the phonon frequency) [17]. The importance of a vari- 
ational phonon basis in predicting accurate results has 
been proved for a two-site system over the entire range of 
e-ph couplings [17]. Our previous study [19] on the 'two- 
site' double exchange model with a single polaron as a 
function of e-ph couplings reveals that the ferromagnetic 
(FM) to antiferromagnetic (AFM) crossover does not co- 
incide always with the polaron crossover. If it occurs for 
some suitable value of J/t, large changes in the effective 
hopping and lattice distortion occur in the crossover re- 
gion on application of the magnetic field. In rcf. [19] a 
possibility of an FM insulating state in between the FM 
large polaron and AFM small polaron state has also been 
observed for small values of J. In this paper we consider 
a many site and many electron double exchange system 
and study the role of the superexchange and e-ph inter- 
action on the magnetic crossover. 

The Hubbard Holstein Hamiltonian in presence of dou- 
ble exchange and superexchange interaction is given by 
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where Ci a (c\ a ) is the annihilation (creation) operator for 
the electron with spin a at site i and is the cor- 
responding number operator, U is on-site Coulomb re- 
pulsion, bi and b\ are the annihilation and creation op- 
erators, respectively, for the phonons corresponding to 
intramolecular vibrations at site i, g denotes the on-site 
e-ph coupling strength. Si and Sj are the core-spins at 
the site i and j respectively, 9 is the angle between the 
core-spins St and Sj and J is the AFM-SE interaction be- 
tween the neighbouring core-spins. The transfer hopping 
integral (t) is modified to tcos(|) by the strong Hund's 
coupling between the spins of the core electron and itiner- 
ant electron [1] . The localized spins are treated classically 
here. It may be mentioned that many workers [6,8] have 
followed the single orbital approximation on the ground 
state as the Jahn- Teller (JT) effect (static and dynamic) 
will split the e g double degeneracy [20] and the mobile 
e g electron would occupy the lower energy orbital at low 
temperature. The present assumption of single orbital is 
expected to be reasonable for doping regime (x <0.4-0.5) 
where CMR occurs. 

For a general value of e-ph coupling, the spread and 
depth of lattice deformations can be studied using Modi- 
fied Lang-Firsov(MLF) transformation. To treat the lat- 
tice deformations that produced at different sites around 
the electron variationally we use the MLF transformation 
[15], 



H = e H He 
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with Ao and Ai represent the lattice deformations at 
the electron site and next-nearest neighbour sites, re- 
spectively. When Ao=5 and Ai=0, MLF transforma- 
tion reduces to Lang-Firsov(LF) transformation [21]. At 
zero temperature the simplest procedure is to make zero 
phonon averaging to obtain an effective polaronic Hamil- 
tonian. But for a many electron system it is better to use 
a two-phonon coherent state [22] for phonon averaging 



) ph = exp 



/ph 



(4) 



a is the squeezing parameter and treated variationally. a 
is nonzero and its effect becomes significant only for in- 
termediate e-ph coupling and finite carrier concentration 
[15]. Increasing a enhances the overlapping of the phonon 
wavefuntions at the nearest-neighbor sites, hence, in- 
creases the polaronic hopping. For a many electron sys- 
tem, consideration of a lowers the ground state energy 
(for intermediate coupling) and smoothens the polaron 
crossover compared to those obtained by MLF and zero 
phonon averaging. 

The MLF-transformcd Hamiltonian, averaged over the 
squeezed phonon state, yields the effective polaronic 
Hamiltonian as 



H cS = £ £p n t,<? ~ t P cos(-) ]T c\ a c 3a 

1,(7 <i.j~>,(J 
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+ V 2 £ n i+s n l+ s, + Nujq sinh 2 (2a) (5) 

where e p is the polaronic self-energy, t p is the polaronic 
hopping, U e g is the effective on-site interaction, V\ and V 2 
are interactions between polarons at nearest-neighbour 
and next nearest-neighbour sites, respectively. These 
quantities are obtained as 

€ p = -(j (2.9 - A ) A + zu \\ 
t p = t exp [exp (-4a){-(A - Ax) 2 - (z - l^ 2 }] 
E/cff = U - 2[w (2.9 - A ) A - zloqXi 2 ] 
Vi = -2(g - A )w Ai 

V 2 = cuoXj (6) 

N is the total number of sites in the system. 

The ground state energy (Eq) of the system per site 
is obtained in the framework of Hartree-Fock approxima- 
tion 



E G , ,0 U cS 2 J 2 

— = e p x e - z p ipcos(-) + — x e + -zS cos9 

+ Vizx 2 e + V 2 z'x 2 e + w sinh 2 (2a) 



(7) 



where z is the number of nearest neighbours, x e is the 
number of electrons per site, 



and p = (ct ffC j, a ) = -j- J2 l^q 



zN 



7q- = J^j e* 1 ?'^ - ^') and Sd (=2) is the spin degeneracy. 
For simplicity we choose a square density of states, then 
one obtains at zero temperature p = ^f(2 — x e ) [15]. The 
polaronic variational parameters Ao, Ai and a as well 
as the magnetic variational parameter 9 are determined 
from the minimization of the ground state energy. For 9 
a simple analytical expression is obtained as, 



9 ptp 



COS 2-2J52 



(for nonzero solution of 9) 



(8) 



For pt p >2JS 2 , 9=0 which corresponds to FM state. In 
this work we will limit to small values of t (t < ujq) so 
that the error encountered in phonon averaging is small 
and does not change our results significantly. We have 
considered lattice deformations only up to the nearest- 
neighbour sites (around the electron) to avoid too many 
variational parameters and because of the fact that for 
t < uiq the inclusion of variational parameters describ- 
ing lattice distortions at distant sites do not change the 
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qualitative behaviour of the results. For numerical cal- 
culation in this paper we consider z and z are 4 and 
12, respectively and energy parameters t, U and J are 
expressed in unit of o>o=l-0. 

In Fig. 1 we plot the angle (9) between the core 
spins on nearest-neighbor sites, the effective hopping (£de 
=£ p cos(0/2)) of the itinerant electron and the on-site lat- 
tice deformation (Xo/g) as a function of e-ph coupling (g) 
for different values of J and £. Figs. 1(a) and 1(b) (for 
i=0.7) show that for small g the ground state is the FM 
with large polaron as carriers. The signature of large po- 
laron is evident from weak polaronic reduction in hopping 
and the reduced value of Xo/g. For large g the ground 
state is the AFM state (9=tt) with small polarons as car- 
riers. The above result may be understood considering 
the fact that increasing e-ph coupling leads to a polaronic 
reduction of the kinetic energy which in turn destroys 
the FM state since the stability of the FM state requires 
pt p > 2JS 2 . Similar behavior has also been obtained by 
recent numerical studies [23]. Fig. 1 shows that the mag- 
netic transition occurs simultaneously with the large to 
small polaron crossover within the present model. With 
increasing J the magnetic transition as well as the po- 
laron crossover shifts towards lower values of e-ph cou- 
pling and becomes smoother. Similar effect is also ob- 
tained by decreasing the value of £ (Figs. 1(c) and 1(d)). 
For £=0.5 and higher values of J (JS 2 — 0.1) the canted 
magnetic state is stable for weak e-ph coupling and the 
transition from the canted state to the AFM state with 
increasing g is very smooth. Fig. 1 also shows that the 
nature of the magnetic transition depends on both t and 
J. 
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FIG. 1. Plot of Xo/g, £de/£ and 9 (in radian) for x e =0.3, 
U=1.0 as a function of g. (a) £=0.7 and JS 2 =0.035; (b) 
£=0.7 and JS 2 =0.07; (c) £=0.5 and JS 2 =0.05 and (d) £=0.5 
and JS 2 =0.1. 
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FIG. 2. (a) Plot of 8, £de/£, Xo/g and Ai/Ao as a function 
of electron concentration (x e ) for £=0.5, £7=1.0, JS' 2 =0.016 
and <7=1.8. Inset figure for higher values of JS 2 and g. (b) 
Plot of 6, £de / 1 vs. x e with £=0.5, JS 2 =0.01 for different 
values of e-ph coupling. 

In Fig. 2(a) the variations of A , A 1; £ DE and 9 with 
electron concentration (x e ) are shown for different values 
of g and J. For low electron concentrations, the ground 
state is AFM and the carriers are small polarons. As the 
electron concentration increases Ao decreases while Ai in- 
creases depicting a small to large polaron crossover. For 
intermediate electron concentrations the ground state is 
FM with large polaronic carriers. The Hamiltonian (1), 
that we considered, has a particlc-holc symmetry. This 
is reflected in Fig. 2(a) where all the physical quantities 
(presented in the figure) arc symmetric with respect to 
half filling (x e =l). In Fig. 2(b) we have plotted 9 and 
tjjE against x e for different values of e-ph coupling. In 
absence of e-ph coupling (g=0) spin canting (0 < 9 < it) 
starts at the smallest carrier concentration and the cant- 
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ing increases with x e until the material becomes ferro- 
magnetic at a low value of critical concentration x cr [24] . 
With increasing g the AFM small polaronic state remains 
the ground state for a range of low carrier concentration 
and the value of x cr increases. For JS 2 /i=0.02, which is 
reasonable for manganites, the AFM-FM transition oc- 
curs around x e ~0.2 for .g=1.8. Fig. 2(b) shows that the 
region of AFM small polaron state extends at the ex- 
pense of the FM large polaronic region with increasing 
J or g (also evident in the inset of Fig. 2(a)). This is 
simply due to the fact that the FM large polaron state is 
destroyed by increasing the AFM interaction as well as 
by the suppression of the kinetic energy with increase of 
e-ph coupling. Hence the combined role of e-ph coupling 
and SE interaction determines the physics of the system. 

It may be mentioned that near half filling (a; e ~l) the 
system is susceptible to show charge ordering depending 
on the choice of the parameters. However, we have not 
addressed this issue here, because our main interest in 
this paper is in the region from low to intermediate filling 
where a transition from an AFM small polaronic state 
to FM large polaron state occurs as the density of the 
carriers increases and the e-ph coupling plays a crucial 
role in this transition. 
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FIG. 3. Bandwidth (W = zt p cos 9/2) as a function of JS 2 
for different values of g with t=0.5, [7=1.0 and x e —0.3. 

In manganites there are SE interactions between near- 
est neighbor e g electrons. This interaction originates 
from the virtual hopping of e g electrons and is propor- 
tional to the square of the hopping matrix element and 
inversely proportional to the increase in energy (<~ on-site 
Coulomb repulsion C/ C ff) m the intermediate state when 
two electrons occupy the same site. The SE interaction 
between the e g electrons determines the orbital order- 
ing temperature in undoped manganite [25]. A relevant 
question may be asked whether this SE interaction will be 
reduced by the small polaron formation as the polaronic 
hopping is suppressed. One of us [26] has shown earlier in 
the context of Hubbard model that the SE interaction is 



not reduced by the polaron formation because in virtual 
hopping the lattice deformations are not transferred from 
site to site. Rather a slight decrease in the value of the 
on-site repulsion due to polaron formation may slightly 
increase the SE interaction between e g electrons. 

In Fig. 3 we plot the electronic effective half-bandwidth 
(W) as a function of the AFM interaction for different 
(intermediate) values of e-ph coupling. It is seen that the 
bandwidth has a unique dependence on J also. In the 
metallic FM region the bandwidth is independent of J; 
whereas in the canted state the bandwidth is suppressed 
drastically with increasing J. Narrowing of the electronic 
bandwidth is considered to be one of the origins of the 
suppression effect of the transition temperature (T c ) for 
doped manganites. Strong e-ph interaction favours the 
band narrowing due to the small polaron formation. The 
AFM-SE interaction may be considered as another can- 
didate for band narrowing mechanism. Our mean-field 
result is qualitatively consistent with the Monte Carlo 
study of FM Kondo lattice model [27]. Fig. 3 also shows 
that increasing e-ph coupling extends the AFM small po- 
laron phase with negligible bandwidth in the JS 2 space, 
similar to that observed as a function of electron concen- 
tration in Fig. 2(a). 

In summary, we have studied the role of SE interaction 
(J) in a many-site DE-Holstcin model within the mean- 
field approximations using the variational phonon basis. 
The polaron crossover and the magnetic transition arc 
studied as a function of e-ph coupling (g) and electron 
concentration for different values of J. The suppression 
of the FM phase and the narrowing of the bandwidth 
with both J and g are observed. The small to large 
polaron crossover as well as the AFM-FM transition is 
found to occur at a higher electronic concentration as J 
or 5 increases. 
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